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Abstract 

Predictions of the extended de Gennes-Fisher local-functional theory for the universal scaling 
functions of interfacial tensions near critical endpoints are compared with experimental data. Vari- 
ous observations of the binary mixture isobutyric acid + water are correlated to facilitate an analysis 
of the experiments of Nagarajan, Webb and Widom who observed the vapor-liquid interfacial ten- 
sion as a function of both temperature and density. Antonow's rule is confirmed and, with the aid 
of previously studied universal amplitude ratios, the crucial analytic "background" contribution to 
the surface tension near the endpoint is estimated. The residual singular behavior thus uncovered 
is consistent with the theoretical scaling predictions and confirms the expected lack of symmetry 
in (T — T c ). A searching test of theory, however, demands more precise and extensive experiments; 
furthermore, the analysis highlights, a previously noted but surprising, three-fold discrepancy in 
the magnitude of the surface tension of isobutyric acid + water relative to other systems. 

Keywords: interfacial tensions; critical endpoints; local-functional theories; near-critical binary fluids; uni- 
versal critical amplitude ratios 
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I. INTRODUCTION 



On passing through a critical endpoint, a binary liquid mixture of, say, species B and C, 
in the presence of the vapor a, exhibits phase separation below T c into two phases j3 and 7 
while above T c there is only one homogeneous liquid phase (3^ [1-4]. The 7 phase, rich in 
C, may be supposed to have the higher density so that in a gravitational field it lies below 
the P phase. In fact there are binary mixtures which display lower consolute points so that 
phase separation occurs above T c and mixing below T c ; but since there are found to be no 
basic differences in criticality, only upper consolute points will be considered explicitly. 

Renormalization group theory has recently shown that the critical behavior realized at a 
critical endpoint lies in the same universality class as on the critical locus that results when 
the pressure is increased so that the vapor phase a is fully suppressed [5, 6]. Nevertheless, 
near a critical endpoint, new bulk and interfacial singularities arise [7, 8]. These are of 
particular interest in the interfacial tensions £g| 7 , £ Q |/3) S a | 7 , and T; a \p 7 corresponding to 
the various interfaces indicated by the subscripts [2-4, 7, 8]. Indeed, except for the first, 
namely E /3 | 7 (T), all these tensions are functions both of the temperature, T, and of the 
overall density, p (or, equally, of the overall composition). 

In terms of the standard notation, with t = (T — T c )/T c , the behavior of the critical 
surface tension E^^T) that vanishes at the critical endpoint can be described by 

£/3| 7 (T) « K\t\», t^O-, (h = 0), (1.1) 

where K is a nonuniversal critical amplitude while in d (< 4) spatial dimensions the expo- 
nent satisfies the scaling relations 

/i = 2/3 + 7 — ^ = 2 — a — v, 

= (3{8 + 1) - v = (d - (1.2) 

first derived, specifically for the surface tension, by Widom [9, 10]. The symbol h in paren- 
theses in (1.1) denotes the ordering field [conjugate to the order parameter ~(p — p c )] that 
is introduced in the scaling description [2, 3, 7, 8, 11, 12]. It is an analytic combination of 
the temperature T and the chemical potentials, \xb and pc, or, e.g., the chemical potential 
difference, Ap = p c ~ Hb, and the pressure, p, (all these being thermodynamic fields in 
the standard sense); furthermore, h is defined so that it vanishes at the critical endpoint, 
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specified, say, by (T c ,p c ,p c ), and is identically zero on the whole 7) coexistence surface 
(or phase boundary) in the (T,[i B -,Hc) thermodynamic space: see, e.g., Fig. 1 of I [12]. 

The other surface tensions between the noncritical vapor or spectator phase a and the 
liquid phases /3, 7, and (3^ share a common "background" contribution, S (T, h), which 
varies analytically and does not vanish at criticality. Thus we can write [7, 8] 

E^CT) » X + |^ + £ o (T,0), t-»0+, (h = 0), 
X a \p(T) « X-|^ + E o (T,0), i->0- (/i = 0-), (1.3) 

where and fT~ are again nonuniversal amplitudes. When, as generally expected [1], 
the wetting temperature, T w , lies below T c , the intermediate j3 phase spreads over the 
interface; then Antonow's rule [1, 12] holds, that is, 

S Q , 7 (T) = E a \ P (T) + S / 3| 7 (T), (h = 0) (1.4) 

which by (1.1) and (1.3) implies 

E ah (T)^(K + K-)\t\v + i: (T,0), i->0-, /i = 0+. (1.5) 

More generally, on the whole thermodynamic surface bounding the a phase, we have 

E(T,h) = £E(T,h) + E (T,h), (1.6) 

where AS(T, /i) and So(7 1 , /i) are the singular and regular (or analytic) parts of the surface 
tension, respectively. Above T c one has E = £ a |/j 7 , while below T c one has £ = £ a | 7 when, 
by convention, h > but £ = when /i < 0. 

According to general scaling principles the singular part can be written asymptotically as 

&E(T,h) w K\t\"S^[m(T,h)], m = M(T,h)/B\t\ 13 , (1.7) 

where the superscripts + and — stand for t < 0, respectively: see I. Because in experiments 
the density is more readily accessible than the chemical potentials the argument m of the 
scaling functions, S^(«), has been chosen proportional to the order parameter which, in 
leading order, may be taken as M = (p — p c ). The coefficient B represents the coexistence 
curve amplitude according to M ~ B\t\P . The nonuniversal amplitudes B and K are 
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introduced in (1.7) to make the scaling functions S^(«) and their arguments dimensionless. 
Then the S^(«) are expected to be universal with, by virtue of Antonow's rule (1.4), 

S M (1)-S M (-1) = 1. (1.8) 

Owing to the technical and conceptual difficulties arising from the fluctuations of capillary 
waves, the renormalization group approach has not so far been successful in calculating the 
S^(m). Instead, advances have been made with the aid of phenomeno logical theories based 
on local-functional concepts. The pioneering studies, both theoretical and experimental, 
have been made by Widom and his co-workers [2-4]. However, as noted in [7, 8], their 
theory predicts a "correction term" varying as |t| 7 which becomes more singular than the 
leading \t\^ term when d > 3 — r\ as is relevant in real three-dimensional systems. This 
defect was remedied in the EdGF or "extended de Gennes-Fisher" theory proposed in [7, 8] 
and implemented recently in I [12] which forms the basis for the present report. By using 
an accurate representation of the equation of state near criticality [11], the EdGF scaling 
functions S^(m) have been calculated explicitly and presented in I. On evaluating the scaling 
functions in zero field, one obtains estimates for the Fisher- Upton ratios [7, 8], namely, 

P = (K + + K~)/K = 0.137 5 ±0.002, Q = K + / K~ = -0.834 ± 0.002. (1.9) 

These ratios have been studied experimentally by Woermann and coworkers [13, 14] for 
the binary fluid 2,6-dimethyl pyridine + water: their results are discussed in Sec. II below. 

To go further and test the basic theoretical predictions for the interfacial scaling functions 
S^(fh) away from the h = axis, the notable experimental data of Nagarajan, Webb, and 
Widom (NWW) [4] for mixtures of isobutyric acid and water offer, to our knowledge, the 
only available opportunity. For the same mixture, Howland, Wong, and Knobler (HWK) [15] 
measured the critical surface tension E ( g| 7 (T) while Greer [16] measured the densities on the 
coexistence curve. These two experiments provide what prove to be valuable consistency 
checks and calibrations for the NWW data as will be seen below. 

It may be noted that there are interesting surface-tension experiments on other quasi- 
binary and binary mixtures, such as n-octadecane and n-nonadecane in ethane [17], 2,5- 
lutidine in water [18], and 2-butoxyethanol in water [19]. However, our attention will focus on 
the isobutyric acid + water system since only in this case are there sufficient measurements of 
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the surface tensions near the critical endpoint to warrant an attempt to extract the interfacial 
scaling functions. 

Since the surface tension vanishes at the critical point as \t\^ with /j, ~ 1.26, and hence 
faster than linearly, the background contribution to the surface tensions proves highly signifi- 
cant as anticipated in I. Thus, in analyzing the experimental data, it is essential to determine 
the background E (T, h) carefully. By virtue of the analyticity of the background, we may 
suppose it can be well represented near the critical endpoint, (t, h) = (0, 0), by an expansion 
of the form 

E (T, h) = E c + E t t + E h h + Z tt t 2 + E th th + Z hh h 2 + ■■■ . (1.10) 

Examination of the NWW experimental data (see their Fig. 8) reveals a definite upward 
curvature in the plots of E vs. T at fixed composition; but the theoretical results for the sin- 
gular part AE(T, 0) alone shows the opposite, downward curvature above T c : see I Fig. 5(a) 
and the negative value of Q in (1.9). To capture this dominant behavior, the expansion (1.10) 
must contain at least quadratic terms in t. 

In the analysis reported by NWW of their own data, the background was assumed to be 
represent able by the form 

S 6 (T,M) ~S c + S t t + S tM tM. (1.11) 

However, this expression is not smooth in t and h, as expected on general grounds, since 
on the critical isotherm it implies [dYi^/di) ~ E t + E^M while M varies as sgn(/i)|/?,| 1 / <5 
with 1/8 ~ 0.21 which is highly singular. Furthermore, the NWW background assumption 
seems to have led to the conclusion that S^{m) above T c was indistinguishable from S^{fh) 
below T c : see NWW Fig. 10. However, when a suitable background of the form (1.10) is 
used in the analysis, one sees that the data support a difference between S^j(fh) and S^{m) 
as, indeed, predicted theoretically in I and expected quite generally: see Sec. V below. 

Units and Exponents 

For convenience of reference we have set out in Table I the exponent values adopted in the 
analyses reported below. They correspond, of course, to those for the three-dimensional 
Ising universality class [20-24]. 
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Since extensive experimental data will be discussed it is useful to specify here the units 
employed. Thus: (i) mass densities will be measured in g/cm 3 and denoted in these units as 
p = p/(g/cm 3 ); (ii) for surface tensions, E, and the corresponding amplitudes, K, etc., units 
of erg/cm 2 will be employed and indicated by K, etc.; (iii) since the reduced free energy 
density will be taken in cm" 3 units, appropriate units for the field h, conjugate to the order 
parameter M = (p — p c ), are g -1 indicated correspondingly by h. The units for various 
critical amplitudes, B, C ± , and the coefficients T, th in (1.10), etc. follow similarly. 

Outline 

The balance of this article is set out as follows. In Sec. II, the theoretical predictions (1.9) for 
the universal ratios P and Q are discussed in the light of the experiments of Woermann and 
coworkers [13, 14]. The density data of Greer [16] along the coexistence curve of isobutyric 
acid + water and the corresponding critical surface tension data of HWK [15] are analyzed 
in Sec. Ill to determine the critical amplitudes B and K. These results are used to calibrate 
the NWW data which are much more extensive but not as precise and, thus, harder to 
analyze with confidence on their own. However, the NWW observations are found to be 
fully consistent with the earlier measurements and, furthermore, confirm the validity of 
Antonow's law in the temperature range studied: see Fig. 2 below. 

On this basis the analysis can be carried forward to determine the crucial background 
term, E (i, h). The results along the coexistence curve and on the critical isochore above T c 
(i.e., for zero field, h — 0) are reported in Sec. IV. A and displayed in Fig. 3. To estimate 
E (t, h) on the critical isotherm, T = T c (t = 0) for nonzero h, it is necessary to re-express 
the NWW observations at constant density in terms of the field. For this purpose one needs 
information regarding the bulk equation of state that goes beyond the value of B. In the 
absence of direct bulk measurement on the critical isotherm, one can progress by appealing 
to bulk universal amplitude ratios (known numerically with appreciable reliability [20-22]), 
combined with the further universal critical ratio S + that relates the surface tension ampli- 
tude K to the correlation length [20]: see (A.l) in the Appendix where the relevant theory 
is reviewed and its application explained. 

By this route one can estimate the amplitude C + for the ordering susceptibility and 
thence, as explained in Sees. IV.B, IV. C, and the Appendix, the necessary conversions can 
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be implemented. Thereby the field-dependence of the background may be gauged. The 
singular part, AH(t,h), then follows by subtraction although, unfortunately, with rather 
limited precision that precludes, for example, any direct determination of the surface tension 
amplitudes, K> and on the critical isotherms. Nevertheless, in Sec. V, the theoretically 
predicted universal scaling functions are compared in suitable plots with the experimental 
surface tension data of NWW. To the extent that consistency is well established, the results 
are encouraging: however, significantly more precise and extensive data close to criticality 
will be required to enable sharper tests of the theory. 

The analysis, furthermore, encounters anew a puzzling dilemma concerning the degree to 
which the dimensionless surface-tension ratio, S + , for the isobutyric acid + water system 
actually conforms to the expected value in light of independent observations of the critical 
scattering and the values found for other systems. The issue is brought up in Sec. IV.B and 
pursued in some detail in the Appendix although without resolution. The considerations 
presented highlight the need for further experiments and, perhaps, for further theoretical 
developments. 

The article is summarized briefly in Sec. VI. 

II. LUTIDINE AND WATER: UNIVERSAL RATIOS P AND Q 

Predictions for the surface tension ratios P and Q defined in (1.9) have been tested ex- 
perimentally by Mainzer-Althof and Woermann (MW) [14]. A binary liquid mixture of 
2,6-dimethyl pyridine [2,6-(CH3)2(C 5 H 3 N)], which is also known as 2,6-lutidine, and water 
was prepared at the critical composition and used in the experiment. The surface ten- 
sions S a i | g(T), E Q | 7 (T), and E a i i g T (T') were measured, for the reduced temperature range 
—0.028 <t< 0.019, by the inverted pendant drop (or rising bubble) method: analysis of the 
contour line of the bubbles gives the surface tension. Also, the mass densities of the liquid 
phases ppiT), p 7 (T), and Ppj(T) were measured. For the critical surface tension, MW re- 
analyzed the experimental data of Kreuser and Woermann [13] using the new MW density 
data. 

In the MW analysis of the surface tension data, the background (1.10) was employed 
truncated at the linear term, S t t. At first MW fitted the data to (1.3) and (1.5) by 
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allowing separate background terms for Y, a \p, E a | 7 , and ^ a \/3-y (i-e., different sets of E c and 
E t ) and varying E c , E t , K ± , and (if + if - ). For the exponent \i they assumed the value 
1.26. The results did not agree with the original Fisher-Upton (FU) values P ~ 0.1 2 and 
Q ~ —0.83 [7, 8]. Furthermore, it appeared that Antonow's rule was not obeyed. However, 
when the analysis was performed with, as demanded theoretically, a common background, 
E (T), Antonow's rule was found to be satisfied (as expected): this procedure yielded 

P cxpt ~ 0.00 ± 0.08, Q expt ~-1.0±0.2. (2.1) 

The P ratio is evidently not quite consistent with the FU value but, significantly, the Q 
value is consistent, both in sign and magnitude. 

Other fitting procedures were also explored by MW by imposing the theoretical value of 
either K + / K [= PQ/(1 + Q)\ or K~ / K [= P/(l + Q)\ and using a common background. 
Averaging the central values of two different MW fits and extending error estimates to fully 
cover the fits yield 

i"expt — O.ll^QQg, Qcxpt — 0.835^q.12 - (2-2) 

These values are consistent both with (2.1) and, now, with the original FU estimates. Fur- 
thermore, the more recent improved EdGF estimates for P and Q from I, quoted above in 
(1.9), agree remarkably well with these albeit biased fits to the experimental data. 

It is appropriate to mention some earlier work. Pegg, Goh, Scott, and Knobler [17] 
studied quasi-binary mixtures of n-octadecane and n-nonadecane in ethane and measured 
the surface tensions through and near both the upper and the lower critical endpoints that 
occur in the vicinity of the hidden tricritical point that arises in this system. Their data, 
although limited, certainly support Antonow's rule and, more importantly for us, suggest 
strongly that the ratio K + / K is negative and of order unity in accordance with the EdGF 
predictions. However, the upper and lower endpoints are less than 0.2 °C apart and, mainly, 
for that reason, the 8 data points for the surface tensions E^| 7 (and the 3 for E a | 7 ) are not 
sufficient to warrant quantitative study. 

The surface tension of the water and 2,5-lutidine system has been measured by Privat 
and co-workers [18] at the critical composition of the lower critical consolute point as well 
as off that critical composition. The authors claim to verify Antonow's rule. However, 
their quoted fits for the amplitudes K, K + , and if - seem problematical relative to the 
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graphical presentation of their data. Thus, and especially in the light of the MW experiments 
on the closely related 2,6-lutidine system, this work cannot be regarded as yielding useful 
experimental estimates for the amplitude ratios P and Q. 

III. ISOBUTYRIC ACID AND WATER: NONUNIVERSAL AMPLITUDES 

A. Coexistence Curve Amplitude 

The coexistence curve of isobutyric acid + water has been studied carefully by Greer [16], 
who measured the coexisting mass densities, pp{T) and p 7 (T), with a precision of 20 ppm 
from 3.5 K below the critical temperature at intervals as small as 5 mK (with a precision 
of ±1 mK) for a sample prepared at close to the critical composition. Two runs were made 
and the data are recorded in Appendix A of Ref. 16. Greer's analysis demonstrated that 
an optimal choice of order parameter (as judged in relation to the number density of a pure 
fluid in a lattice-gas representation) was the volume fraction of one component, say 4>b- 
The surface tension experiments of NWW, which are of principal concern to us, however, 
also used the mass density in the vicinity of the critical endpoint as the primary controlled 
observable. Accordingly, we have reanalyzed Greer's data using the mass difference, p — p c , 
as the order parameter. This choice, of course, affects the real (and apparent) magnitudes 
of the various correction terms; but these will play no more than an auxiliary role in our 
analysis. Further discussion of fits for the volume fraction difference, A0, are given in [25]: 
one learns that a range of reasonable fits to the data serve to determine the critical point to 
within ±10 to 15 mK while the leading critical amplitude B can be found with a precision 
of no better than ±1 %. 

More specifically, we have fitted the data to 

±Ap = \ [p 7 (T) - p p (T)] = B\tf [1 + b e \t\ e + M], (3.1) 

with t = (T — T c )/T c as usual, but, as a useful crosscheck, also in terms of the asymptotically 
equivalent variable, 

t' = l-T c /T = t/(l + t), (3.2) 
as in a previous study [26]. Using primes to denote the amplitudes fitted with t' we find, 
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from the first run [25] , 

T c ~ 25.996 7 °C, B ~ B' ~ 0.03104, 

b e ~ 6' e ~ -1.73, 6i ~ -9.783, 6; ~ -9.455. (3.3) 

The difference b[ — b\ ~ 0.328 is close to /3 as it should be for consistency in light of (3.2). 

The exponent 9 ~ 0.5 (See Table I) in (3.1) specifies the leading correction to scaling. 
General scaling considerations [27, 28] indicate that a further correction term, |t| 2/3 , should 
also be present which, since 2(3 ~ 0.65 < 1, should in fact dominate the linear term. The 
corresponding amplitude 62/3 may well be rather small: nevertheless, the data cannot resolve 
two such terms so that the fitted values of the coefficients b g and b' e must be regarded as no 
more than effective amplitudes. 

On setting bg = b\ = in (3.1) Greer found B = 0.0265 ± 0.0005; however, on allowing 
for b e ^ the value 0.0315 ± 0.0025 resulted; this is fully consistent with (3.3). Repeating 
the same procedure for the second experimental run yields 

T c ~ 25.969 6 °C, B ~ S' ~ 0.03175, 

b e ~ b' e ~ -2.635, 61 ~ -14.563, b[ ~ -14.234. (3.4) 

The difference b[ — bi~ 0.329 is again close to (3. Clearly, the deviations of the fitted values 
here from those in (3.3) provide a measure of the accuracy available in estimating T c and B. 

B. Comparison with NWW Coexistence Curve Data 

We may now use the fits to (3.1) to calibrate the experiments on the same system by NWW 
who prepared mixtures of isobutyric acid + water at various compositions and determined 
the surface tensions by optically measuring the wavelength of surface waves generated by 
a transducer. In Table 1 of NWW, about 140 (T, p, S) data points are presented: these 
are what we study here. We have extrapolated some of the observations to the coexistence 
curve; but, as indicated briefly in Sec. IV.A, only very small changes in the values of p and 
E are entailed. 

To proceed, we first consider the diameter of the coexistence curve in the density variable, 
namely, 

p(T) = i( P/3 + p 7 ), (3.5) 
10 



since, as mentioned, the density was directly observed by NWW. Owing to the mixing of 
thermodynamic fields near criticality in fluids, one expects a dominant singular term |t| 1_a 
to appear in the diameter [27-31]. However, the NWW data do not reveal a signature 
of any such singular term. Rather, the coexistence curve appears almost symmetric in the 
(p, T) plane [32]. Thus the diameter may be well represented by a constant as 



As discussed by Greer [16] , the critical temperature of the isobutyric acid + water system 
is particularly sensitive to ionic impurities [33] and, indeed, experience shows that the ob- 
served critical temperatures of nominally the same binary fluid mixtures typically vary from 
experiment to experiment by amounts exceeding the stated errors. For fitting the NWW 
data, therefore, we have adopted their value of T c = 26.310 ± 0.001 °C, even though it lies 
outside the uncertainty limits implied by (3.3) and (3.4). As seen in Fig. 1, by using (3.1) 
and (3.6) with the amplitude estimate 



the NWW data can be represented rather satisfactorily. This demonstrates the consistency 
of their observations with other careful studies and establishes a reliable estimate for the 
amplitude B which will be employed below. 

It may be remarked that NWW advocated the use in theoretical analysis of a (con- 
ventionally defined) volume fraction, 0, in place of the density, p [34]. Study of their data 
reveals, however, that the dependence of on p is both surprisingly nonlinear and signif- 
icantly temperature dependent. Furthermore, contrary to the expectations of NWW (see 
p. 5781) and Greer, the coexistence curve appears to be more symmetric and regularly be- 
haved in terms of p than of 0. Accordingly, we have accepted the directly observed variable 
p as the order parameter in all the following analysis. 

C. Critical Surface- Tension Amplitude K 

In order to obtain the critical amplitude K as defined for the surface tension in (1.1), the 
capillary-rise data obtained by HWK [15] for isobutyric acid + water will be employed. In 
their Fig. 8, one can clearly see one point that deviates significantly from the general trends, 



f{T) ~ p c ~ 0.9936. 



(3.6) 



B = 0.0314 ±0.0004, 



(3.7) 
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namely, that corresponding in HWK Table VII to T c — T = 1.088 K or t ~ -0.0036, and 
= 0.0191. This point has thus been omitted in the present analysis. 
Since only eight data points are available, the correction-to-scaling terms are not expected 
to play a detectable role. First, therefore, the data have been fitted to the form 

Ep\ y = E c + K\t\' i , (3.8) 

where, of course, S c = is anticipated. The reported value of T c is 299 K [15]; but we 
have not found in the account of HWK a satisfying description as to how this value was 
determined. We therefore adopted 

T c = 299.00 K = 25.85 °C, (3.9) 

provisionally but have also examined how estimates for K are affected by changes in the 
assigned value of T c . 

A weighted least-squares fit then yields 

t c = -0.002 ±0.008, K = 34.6 ±5.7, (3.10) 

where the uncertainties quoted here and below are twice the standard deviation. By imposing 
S c = 0, which is clearly consistent with the data, one finds 

K = 32.8 ±1.7, (3.11) 

where the central estimate has shifted by 5% while the uncertainty is reduced by 30% in 
comparison with (3.10). Repeating the analysis with S c = but using t' — t/(l + 1) as a 
variable [26], we find K = 32.6 ± 1.7 in good agreement. 

As a further check one can calculate E^/l^ from each data point. This ratio approaches 
K when t — > 0, and an average yields K ~ 32.7, suggesting that the corrections to the pure 
power law are relatively small. 

To gauge the sensitivity to the imposed value of T c , we changed the value by ±0.1 K. 
Considering that the temperature was controlled to within 0.01 K [15], this variation is fairly 
large. The results (obtained by fixing S c = 0) are K = 30.6 ± 1.7 for T c = 25.95 °C and 
K = 35.3 ± 1.8 for T c = 25.75 °C. Despite the changes of ±7 %, the overlap of the error 
bars allows consistency with the previous estimates for K. 
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Finally, the data were fitted to E = K | (T/T c ) — 1^ treating both T c and K as parameters. 
Minimizing X 2 yielded K ~ 33.8 5 and T c ~ 25.78 °C. This value for T c is 0.07 K lower 
than the asserted value (3.9) but the amplitude estimate again lies in the indicated range. 
Repeating the procedure with E = K' |1 — (T C /T)| M , yields the same value for T c but K' ~ 
33.5, about 1% smaller. Overall we believe that 

K = 33.7 ±2.5, (3.12) 

summarizes the evidence conservatively [35], although a more precise knowledge of T c in the 
HWK experiments could be valuable. 

D. Comparison with NWW Critical Surface- Tension Data 

It is informative to compare the experiments performed by HWK with those of NWW to 
check, in particular, for mutual consistency: see Fig. 2. For convenience NWW data points 
on the coexistence curve have been numbered consecutively from 1 to 15. Except for point 8, 
other points joined by parallelograms were measured at the same temperature. Points 1-7 
lie on the high density side of the coexistence curve and hence represent E a | 7 (T, h — 0+), 
while points 9-15 correspondingly lie on the low density side and represent E^^T, h — 0— ). 

To explain the significance of the fixed-T parallelograms, consider, as an example, points 2 
and 14. By Antonow's rule we have E a | 7 = Eq,^ + E^ and so by using E /3 | 7 (T) « K\t\^ 
with the estimate (3.12) and by varying Eq,^ within the uncertainty limits, one can find a 
range of E a | 7 values that are consistent with Antonow's rule. Evidently, the central value 
of the data point 14 corresponds to the lower limit of the data point 2, while the upper 
limit corresponds to the central value of data point 2. These mutually consistent limiting 
pairs of the data points 2 and 14 have been joined by parallel lines. In other words, the 
vertical sides of the parallelogram connecting data points 2 and 14 represent error limits on 
the NWW data that are consistent with the HWK experiment and with Antonow's rule. 
Thus we notice that data point 3 represents a high estimate of Y, a \p while data point 13 is 
relatively low. Despite this (3,13) worst case, it can be concluded that the two experiments 
are mutually consistent and that the estimate K in (3.12) is reliable [35] . Note that it would 
have been very difficult to obtain a reliable value of K from the NWW data alone, owing 
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to the necessity of extrapolating and differencing their (not highly precise) data to obtain 
£g| 7 via Antonow's rule. 



E. Other Surface Tension Amplitudes 

Having determined the two amplitudes K and B, which serve as metrical factors in the 
scaling formulation of the critical-endpoint behavior of the surface tension [see (1.7)], we are 
prepared to take the next step towards analyzing the NWW data. The first issue must 
be to establish the background contribution, S (T, h), by estimating the coefficients E c , S t , 
etc., in (1.10). To that end it proves helpful to have at hand values for the surface tension 
amplitudes K + and K~ as defined in (1.3). Ideally, of course, these should be unambiguously 
determined by the data themselves: but, as already seen in Sec. II, it is unrealistic at this 
stage in the development of experimental techniques to expect to do more than verify- 
optimistically at a fully convincing level — consistency of the theory with the observational 
data. Accordingly, we report here the values 

K + = -23.3 ± 1.8, K~ = 27.9 ±2.1, (3.13) 

that follow from (3.12) on accepting the EdGF calculations for the universal ratios P and Q 
given in (1.9). 

To determine the field-dependent coefficients T, h and E h/l in the background expansion 
(1.10) we will need to examine the data for T,(t,h) on the critical isotherm. The expected 
behavior is 

E c (h) ^Ki\hf A + ^ c + ^ h h + ^ hh h\ (T = T C ) (3.14) 

where, again, values for the amplitudes K> and K<, for h < 0, respectively, will prove of 
importance. To estimate these one may appeal to the EdGF theoretical values reported in 
I (6.9) for the universal ratios 



In the second factor in this expression, which enters because K> and K have different di- 

< 

mensions, the nonuniversal amplitude C + specifies the magnitude above T c of the divergence 
of the basic ordering susceptibility as l/f. Note that the notation for amplitudes used here 
and below accords with that set out in the Appendix of 
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The appropriate value of C + , and of other nonuniversal amplitudes for isobutyric acid + 
water, is taken up in Sec. IV.B when the field-dependence is considered explicitly. At that 
stage, we will make contact with Moldover's extensive analysis [36, 37] of surface tension, 
light scattering, and bulk thermodynamic experiments for a range of single-component and 
binary fluid systems designed to test the hyperscaling (or "two-scale factor universality") 
hypothesis for criticality. Various details and some questions they raise are presented in the 
Appendix. 

IV. DETERMINATION OF THE SURFACE TENSION BACKGROUND 

As discussed in Sec. I, to analyze the NWW surface tension measurements appropriately, 
the background S (T, h) must at least contain quadratic terms in t; it then seems reasonable 
to include quadratic terms in h also as in (1.10). In order to determine the background 
coefficients £ c , S t , etc., we adopt the following strategy. First, the data on the h = axis, 
i.e., for the mixture at the critical composition, will be fitted to obtain E c , S t , and S tt . From 
the data on the critical isotherm, t = 0, the coefficients and T*hh can then be examined. 
Finally, to obtain the cross-term coefficient Ti th , the derivative (<9£/<9t)h may be studied on 
the critical isotherm. 

A. Along the Coexistence Curve 

In zero-field, the surface tension above T c is, following the Introduction, expected to behave 
as 

Eal^T) » K + \t\» + S c + E t f + 5V 2 , (T > T c ). (4.1) 

Below T c , there are two distinct vapor-liquid interfaces and two surface tensions, H a \p and 
S a | 7 , corresponding to the coexisting phases (5 and 7. By Antonow's rule (which as shown 
is satisfied by the NWW data) the difference (£ a | 7 — ^ a \p) is the critical surface tension 
£g| 7 and thus it suffices to study only Yi a \p(T) below T c . The identical expression (4.1) then 
applies except that K + is to be replaced by K~ . 

To proceed, one may first fit ^^^(T) using the data above T c at the very close-to-critical 
composition tabulated by NWW, by regarding K + , E c , S t , and S tt as free parameters. 
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These fits suggest K + > which, in accord with the generally observed upwards curvature 
(compare with I Fig. 5), has the opposite sign to the prediction for K + in (3.13). We 
conclude that the NWW data above T c are not, on their own, of adequate precision to 
reliably distinguish the amplitudes of three terms varying as t, t M , and t 2 ; hence, we cannot 
independently estimate K + . 

However, as indicated in the MW analysis [14] discussed in Sec. II, the situation may be 
improved by imposing the theoretical value of the ratio K + /K; but this is then equivalent 
to imposing K + because, as seen in Sec. Ill, the critical amplitude K has been already 
estimated reliably by comparison with the HWK data. Since E a | y 9 7 (T) and Y, a \p(T) share 
the common background S (T, 0), we likewise impose K~ as given in (3.13). Of course, by 
adopting this procedure we will be primarily checking the consistency of our theory with the 
NWW data which, unfortunately, are not adequate for estimating K + and K~ directly as 
might originally have been hoped. 

Accordingly, by subtracting the singular term in (4.1) and likewise below T c , we can 
generate experimental data for the zero-field background alone: see Fig. 3. Again we use 
the NWW data taken at the (close-to) critical composition above T c and the (close-to) 
low-density side of the coexistence curve. A least-squares fit then yields 

t c ~ 26.284, t t ~ -14.6, t tt ~ 381. (4.2) 

The corresponding plot is represented by the dashed curve in Fig. 3 and seen to be very good. 
This demonstrates that these coefficient estimates are most reasonable although, as will be 
discussed, they prove to be less than fully satisfactory in the near critical region below T c . 

To examine the fit further consider the plots shown in Fig. 2. The dashed curve represents 
the surface tensions H a \p(T) and £ Q | 7 (T) on the coexistence curve as reproduced from the 
amplitudes K and K~ and the background coefficients (4.2): to translate t oc T c — T into the 
density on the coexistence curve, the inner fit described in the caption of Fig. 1 is used. The 
agreement with the central values of the NWW data on the low-density side is encouraging; 
however, the fit reveals significant, systematic discrepancies on the high density side where 
the plot lies outside the observational range for all the points 1-7. 

To complement the overall or "full-range" approach just described and to understand 
the discrepancies uncovered on the high-density side of the coexistence curve, let us re- 
turn to Fig. 2 and focus on the data closest to criticality. In particular, as explained in 
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the construction of the parallelograms in this figure, consistency with the HWK data for 
E ( 3| 7 (T'), together with Antonow's rule [38] requires that the allowed observational uncer- 
tainties quoted for the NWW data should be reduced to correspond to the vertical sides of 
each parallelogram. Accordingly, for T < T c consider the truncated forms 

E^T) ~ X-|^ + E c + E t *, (4.3) 
E a , 7 (T) ~ (K~ + K)\t\» + E c + £ t i, (4.4) 

in which the background has been restricted to a linear term. By trial and error, one can 
discover those values of E c and E 4 that give acceptable fits to those data points lying closest 
to the critical point, specifically, for the pairs (7,9), (6,10), (5,11), (4,12), and (3,13). Slight 
shifts in the nominal value of p c from (3.6) may also be examined but are found to be 
detrimental for changes exceeding Ap c = 10~ 4 . 

The predominant conclusion is that E c , the common value of the surface tensions H a \p, 
E Q i 7 , and ^ a \/3y a ^ criticality, cannot actually be as low as suggested by the full-range fit 
(4.2). The acceptable "critical fits" yield, instead, the range 

E c = 26.294 ±0.004. (4.5) 

Indeed, NWW actually quote the estimate E c ~ 26.30, which is much closer to this "critical 
value" than to (4.2). The corresponding preferred value for the background slope, is found 
to be E t = —17.5; but reasonable fits could be obtained in the range —16.5 > S t > —18.5. 
In all these fits, however, the fitting curve in Fig. 2 became too broad at the larger values of 
\p — p c |, failing to fit the data-point pairs (1,15) and (2,14). This last defect, however, was 
overcome by including in (4.3) and (4.4) the quadratic term E tt t 2 . As the assigned value E tt 
is increased, keeping E c fixed as proves essential, the optimal value for E t correspondingly 
falls. Thus E t = —16.0 and E tt = 250 give very good fits in Fig. 2. However, these particular 
values do not describe well the highest data points above T c in Fig. 3 (for t > 0.015). One 
could well disregard this range since the data extend below T c only down to \t\ ~ 0.011. 
However, one can find values that prove rather satisfactory in fitting all the h — data, i.e., 
on the coexistence curve and at the critical composition above T c , namely, 

E c ~ 26.294, E t ~ -15.0, t tt ~ 400. (4.6) 
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These coefficients yield the solid curve shown in Fig. 3: evidently the fit displays a systematic 
displacement above the central values of the data points but by rather small amounts that 
remain within all the uncertainty ranges. Furthermore, the agreement with the coexistence 
data is now much improved as seen by the solid curve in Fig. 2. Only very minor departures 
from the acceptable uncertainty ranges now arise. Accordingly we will retain the assignment 
(4.6) in all the following analysis. 

At this stage we may conclude that the NWW data are, first, fully consistent with the 
best previous observations of the coexistence curve and the critical surface tension and, 
second, are quite consistent with the EdGF based estimates of the ratios K + / K and K~ / K 
even though, unfortunately, the data cannot critically test these ratios. 

B. On the Critical Isotherm 

It transpires, in retrospect, that NWW did not make observations of the surface tensions 
very close to the critical isotherm T = T C = 26.310 °C. We may, however, presume that the 
isotherms just above and below T c , namely, at T = 26.49 °C and T = 26.05 °C constitute 
rough lower and upper bounds, respectively, on £ c (/i). On the critical isotherm, the surface 
tension should behave as in (3.14) with /i/A ~ 0.806. The optimal assignment of the critical 
point value S c has already been determined in the previous subsection. One might wish 
to determine the amplitudes and K C K directly from the data; but, regrettably there is, 
once again, no hope of that! Rather, with the limited aim of testing only consistency with 
the theory, we will adopt theoretical estimates for these two amplitudes. At this point one 
encounters a rather sharp and surprising dilemma: as explained further in the Appendix, 
two significantly different sets of EdGF-based predictions arise from (3.15) and I (6.9). The 
first set, derives entirely from the two metrical factors K and B, that we have determined for 
isobutyric acid + water, combined with well-confirmed values for various universal critical 
amplitude ratios [20]: as explained in the Appendix this yields 

A: K.%, = (5.0±0.6)xl0" 18 , K% = -(1.5 5 ± 0.1 7 ) x 10~ 18 . (4.7) 

The second set also utilizes the estimate (3.12) for K in (3.15) together with (3.7) for B; 
but instead of using the purely theoretically derived value, namely, 

A: C + = (3.04 ± 0.30) x 10~ 26 , (4.8) 
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for the amplitude C + in (3.15), it employs the much larger estimate 

B: C + = (14.2 ± 0.6) x 10" 26 , (4.9) 

that follows, as explained in the Appendix, from independent light scattering and thermo- 
dynamic experimental data for isobutyric acid + water [36, 37]. This then leads to the 
alternative estimates 

B: Kl = (17.5 ± 1.5) xlO" 18 , K% = -(5.4 ± 0.5) x 1(T 18 , (4.10) 

that, likewise, are much larger than in (4.7). 

Granted the values of K> in the critical point background (3.14), one is left with only the 

< 

coefficients Yu h and T, hh to be determined by fitting. However, since the NWW measurements 
were made in terms of the density, it is necessary to convert the (t, M = p—p c ) data to (t, h) 
using an equation of state. To this end, it is natural to employ the scaled equation of state 
derived in [11] on the basis of the extended sine model. The metric factor B given in (3.7) 
and the appropriate values of C + [in (4.8) or (4.9)] may be used. The data spread over 
the range —4.8 x 10 21 < h < 2.4 x 10 20 [using (4.8)] and, hence, reach more widely on the 
negative h (or low density) side; but for fitting we will focus on the available symmetric 
range. 

On the basis of the A values above we thus find that the assignment 

t h ~ 3.72 xlO" 23 and t hh ~ 6.65 x 10~ 43 , (4.11) 

provides a fit that lies well between the upper and lower limits for \h\ not too large: see 
Fig. 4. Comparable fits can be found using the B values but with, of course, substantially 
different values of and and a different accessible range of h. For the balance of this 
article, however, we will stay with the 'purely theoretical' A estimates (4.7) and (4.8). 

C. Cross Term 

In order to estimate the cross coefficient Ey, in the background (1.10), one should examine 
the derivative (dH/dt)h at T = T c . Since the measurements were made at constant M (or 
density), it is appropriate to invoke the relation 




(4.12) 
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However, owing to the measurement uncertainties, the estimates of (dH/dt)M obtainable 
from the NWW data are extremely noisy. Hence, one cannot seriously attempt to estimate 
T, th . One may reasonably guess, nevertheless, that the geometric mean of E tt and T, hh , 
namely, using (4.6) and (4.11), 

(£«£m) 1/2 = 1.6x10- 20 , (4.13) 

will provide some indication of the possible order of magnitude of the cross-coefficient. How- 
ever, we comment below on the sensitivity of the final comparisons with experiment to the 
value of E t fc. 

V. SCALING FUNCTIONS FOR THE SURFACE TENSION 

Having found, to the somewhat limited degree feasible, a reasonable description of the im- 
portant surface-tension background, we may examine the experimental data for consistency 
with the predicted universal scaling functions S^(m) introduced in (1.7). 

In order to subtract the background E (t, h), from the full surface tension data, it is 
necessary, as before, to convert the variables (t,M — p — p c ) to (t, h) because the exper- 
iments employed the density as a primary variable. Adopting the NWW critical tem- 
perature T c = 26.310 °C and using B from (3.7), the first step is to compute m = 
(p — Pc)/B\t\P. The universal, scaled equation of state m = Q±(h) [for t < 0] may be 
used in the form of the extended sine model of [11] to implement the conversion to h. The 
corresponding h follows via h = Bh\t\ A /C + where, for the present discussion we use (4.8) for 
C + as in Fig. 4. Then, using the background with the coefficients given in (4.6) and (4.11), 
an estimate for the singular part of the surface tension, AE(T, h), can be extracted from the 
NWW data. 

Note that in the first instance, the cross term, H t hth, may simply be dropped. The 
resulting scaling contribution, AE/i^|t| M , is plotted versus m, the scaled density deviation, 
in Fig. 5. 

Since the field-dependence of the background E (t, h) has been estimated only from data 
that lie on the two isotherms closest to criticality and fall within the restricted range \h\ < 
2.4 x 10 20 , we have distinguished, by solid dots versus crosses and plusses, between these 
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observations, corresponding to the restricted range of h, and (selected) further data falling 
outside this range. 

Because the background contribution constitutes such a large part of the total surface 
tension, it is not surprising that the subtraction needed to isolate the singular piece, A£(t, h), 
is subject to relatively large uncertainties. This appears to be the main reason why many of 
the crosses and plusses, derived from the low-density data (beyond the fitted range of h) , fall 
significantly below the scaling loci. For the data inside the range (the solid dots in Fig. 4), 
the agreement between the theoretical and experimental results is moderately encouraging 
considering the uncertainties of the original data and the inherent difficulties of the analysis. 
Certainly, the data of NWW support the expected difference between S^(m) and S^-(m). 

It may be noted, nevertheless, that just as pointed out by NWW, the logarithmic singu- 
larity associated with the complete-wetting transition cannot be detected in these or other 
surface-tension plots. However, this is due not only to the limited precision of the experi- 
mental data but also to the fact that in the usual representations of the data the singularity 
is not really visible even in theoretical plots where it is known analytically to be present: 
compare with Figs. 5 and 6 of I. 

If one also includes a positive cross-coefficient Hth in the background with the magnitude 
given in (4.13), one finds slight changes of less than 1% in the data points plotted in Fig. 5. 
However, when the cross-term is inserted with a negative sign, the majority of the data 
points display changes of ~5% and a few points prove still more sensitive. This situation is 
hardly satisfying but, in light of the experimental challenges and the resulting precision of 
the NWW data, it is clear that our analysis can proceed little further. 

Finally, as pointed out theoretically by Ramos-Gomez and Widom [3] and mentioned 
briefly in I, the surface tension isotherms when examined vs density (oc M) above T c should, 
in general, exhibit crossings. Although hints of such crossings have appeared in some binary 
fluid systems [39-41], they were not seen in the NWW experiments. The primary reason 
seems to be that for the range of densities p > p c explored, the isotherms chosen for ob- 
servation above T c were somewhat too widely spaced. Indeed, as illustrated in Fig. 6, the 
anticipated crossings for isobutyric acid + water can be exhibited on the basis of our fits to 
the NWW data. All that is needed is to choose temperatures above T c spaced apart by 
intervals of no more than 1 °C and to extend the observations to densities slightly further 
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above p c . The crossings, while sensitive to the presence of the background term, should be 
visible even in the singular part of the surface tension, AE(T, p). Note, however, that since 
our fits were based only on the restricted (symmetric) range of h, the plots lying above the 
short horizontal lines in Fig. 6 may be less reliable than at lower densities. 

VI. SUMMARY 

In pioneering experiments to study critical endpoint behavior, Nagarajan, Webb, and Widom 
[4] made extensive measurements of the vapor-liquid interfacial tension of isobutyric acid and 
water mixtures as a function of temperature at various compositions. Since the precision of 
the surface tension data was somewhat limited, we used the observations of the critical surface 
tension made by Howland et al. [15] and the coexistence curve measurements of Greer [16] to 
cross-calibrate the NWW data. We found full consistency between the various observations 
and verified agreement with Antonow's rule near the critical endpoint. In this way reliable 
estimates of the coexistence curve amplitude B and the surface tension amplitude K were 
obtained: see (3.7) and (3.12). Using known universal amplitude ratios [20] and these two 
amplitudes, which can serve as the only required metric factors in a hyperscaling formulation, 
we estimated the nonuniversal susceptibility amplitude C + . Thereby the experimental data 
as a function of density could be transformed in order to uncover the field-dependence of the 
crucial surface-tension background contribution in the form (1.10). Owing to the considerable 
uncertainties, however, we could not resolve the cross-coefficient Sj/j. Nevertheless, even 
without this term, the scaling plots for the singular part of the interfacial tension as derived 
from the experimental data proved consistent with the theoretical predictions based on the 
EdGF theory of I [12]: see Fig. 5. 

On the other hand, as discussed in the Appendix, significant discrepancies come to light 
when the Howland et al. data are compared, via hyperscaling expectations, with evidence 
from other critical systems. Owing to conflicting experimental reports [35, 36], however, it 
seems impossible to resolve this issue and assess its significance without further experiments 
and, perhaps, new theoretical insights. 

More rewarding from a theoretical perspective and less open to questions of precision, the 
critical surface tension measurements of Mainzer-Althof and Woermann [14] on 2,6-dimethyl 
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pyridine + water at the critical composition, show rather satisfactory agreement between the 
experimentally determined universal ratios P and Q [see (1.9)] and the EdGF predictions. 
However, as regards the full scaling aspects of the theory, more detailed, more precise and 
more reliable data are sorely needed to execute more searching tests! 
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APPENDIX: NONUNIVERSAL CRITICAL AMPLITUDES 

In Sec. Ill the coexistence curve amplitude B for isobutyric acid + water (IBA:W) was de- 
termined from the data of Greer [16] and checked against the NWW data yielding (3.7). 
Likewise, analysis of the HWK capillary-rise surface tension data yielded the critical ampli- 
tude K reported in (3.12) that also accords well with the NWW data using forced capillary 
waves: indeed, no larger deviation in (3.12), than, say, ±6.0 (or ±18%) could be reconciled 
with the NWW observations [42]. We stress this point because Moldover [36], in comparing 
these IBAW data with results for many other fluid systems, uncovered an inconsistency in 
the value of the expected-to-be universal amplitude ratio 

S + = K(i+f = 0-377 ± 0.011 with K = K/k B T c . (A.l) 

Here, in the notation of the Appendix and Table I of [20], the amplitude £f specifies the 
magnitude of the divergence of the second- moment correlation length above T c via £i(T) 
£,t /t v ■> while the value quoted for S + is that reported in [20]. 

If we now accept this theoretical result we may proceed to predict the value of for 
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IBA:W. Thence, using the universal ratio [20] 

S = C + K/B 2 £+ = 1.17 ± 0.06, (A.2) 

together with the estimates for K and B, we can compute the susceptibility amplitude C + 
that is needed in Sec. V to establish the h scale in Fig. 4. This route, indeed, yields the A 
value for C + reported in (4.8) and, with the further aid of (3.15) and I (6.9), the values for 
K> and K< given in (4.7). These results, in turn, are incorporated in the plots in Fig. 5 
that serve to test the scaling hypothesis. 

o 

The correlation length amplitude predicted in this fashion is £^ pre = 2.10 ±0.08 A [where 
the HWK value (3.9) has been adopted for T c in (A.l)]. However, as reported by Moldover 
[36], in light scattering experiments by Chu and co-workers [43, 44], the second-moment 
correlation length above T c was evaluated for IBA:W. The fits presented led them to conclude 

o 

v ~ 0.613 and = 3.57 ± 0.07 A. This correlation amplitude exceeds the predicted value, 
£i" P re' by 60% or more! On re-analyzing their data with the higher value v = 0.632 imposed, 

o 

we obtain = 3.09 ± 0.05 A. This is somewhat closer to £i" pi . e although the uncertainty 
limits fail to overlap by 36%. From turbidity measurements, Beysens et al. [45] concluded 

o 

= 3.62 5 ± 0.06 5 A and v ~ 0.630; this value of agrees with Chu et al. and is some 17% 
higher than that from our reanalysis. 

At a basic theoretical level the discrepancies uncovered here are of serious concern. If, for 

o 

concreteness, one accepts £^ cxpt ~ 3.6 ± 0.07 A as the observed value, in accord with these 
experiments, the universal ratio S + in (A.l) would yield a prediction for K of only 12.0, about 
36% of that observed! Conversely, one might be tempted to assert that the experimentally 
determined value of the (supposedly universal) ratio for IBA:W is S^ pt = 1.06±0.03, in place 
of (A.l), the uncertainty indicated here taking into account various sources of experimental 
and theoretical imprecision. 

Similarly, the B values for C + and K% reported in (4.9) and (4.10) follow from the 

< 

observed value of B and the experimental value, £^ cxpt , by using [20] 

Q c = (£+fB 2 /C + = 0.323 6 ± 0.006. (A.3) 

In facing this discrepant situation let us focus first on the theoretical value (A.l) for 
S + (which, indeed, was incorporated into the numerical scaling calculations reported in 
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I and tested here). One might initially notice that 0.377 is significantly higher than the 
various theoretical, RG and Monte Carlo estimates, namely, ~0.21, 0.24, and 0.28, available 
to Moldover in 1985 [36]. In fact (A.l) rests on the analysis [26] of difficult but carefully 
performed Monte Carlo simulations for the nearest-neighbor (d = 3)-dimensional Ising model 
by Hasenbush and Pinn [46] (combined with the well established universal ratio = 
1.96 ± 0.01 [26, 47]). The simulations [46] could approach criticality no closer than \t\ > 
0.015. By contrast, the experimental data for IBA:W lie in the region \t\ = 0.0008 to 0.010. 
Although it seems unlikely, it is conceivable that exceptionally strong correction-to-scaling 
terms in the surface tension of the simple cubic Ising model could dominate very close to T c 
so that (A.l) represents a significant under-estimate of S + . Alternatively, finite-size effects 
in the simulations may be playing a far larger role than appreciated although no evidence to 
suggest this was observed and the issue was by no means neglected [46, 48]. 

On the other hand, it might be that for reasons associated with the lattice structure and 
capillary-wave fluctuations, the standard simple cubic lattice gas/Ising model is inadequate 
for certain real fluid mixtures. Most certainly, the conventional near-neighbor lattice gases 
neglect the long-range power-law van der Waals interactions that lead to interphase potentials 
decaying, normal to an interface, only as 1/z 3 . And while these van der Waals forces play 
only a subdominant role in bulk critical behavior, they are known to be especially significant 
in surface and wetting phenomena [49]. It is, thus, possible that they seriously distort the 
near critical interfacial tension or even destroy the expected universal character of the ratio 
S+. 

At present, however, effective theoretical tools for addressing this interesting issue are not 
apparent. Furthermore, it is not unfair to remark that these various rather deep theoretical 
issues do not directly affect our basic scaling analysis of the NWW observations of surface 
tensions near a critical endpoint. Indeed, in the first instance only the experimentally deter- 
mined critical amplitudes K and B enter: universal or nonuniversal relations between these 
and other critical amplitude do not play a definitive role in the primary scaling issues. 

Finally, we must return to the experimental evidence. It is, indeed, striking that the 
value for S + in (A.l) accords remarkably well with the experimental values assembled by 
Moldover [36] for Ar, Xe, N 2 , O2, CO2, CH 4 and SF 6 and for the binary fluids triethylamine 
+ water, nitroethane + 3-methylpentane, cyclohexane + aniline, methanol + cyclohexane, 
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and even for various polystyrene + C 7 H 14 solutions: see Fig. 1, Table I and also the Note 
added regarding Ref. 86 in [36] and Table I of [37]. Indeed, by and large, Moldover's analysis 
confirms rather generally the hyperscaling hypothesis, both for the interfacial and for the 
bulk thermodynamic properties. (See also in [25].) 

These observations tend to reinforce confidence in the estimate (A.l) for the universal 
surface tension amplitude — at least as regards 'simple' pure fluids and binary mixtures. But 
does the weak acid IBA:W constitute such a 'simple' system? On the one hand, as an 
electrolyte, a population of positive and negative ions will be present in all the liquid phases. 
Furthermore, it is known that the critical point of IBA:W is extremely sensitive to small 
impurities and, especially, to ionic impurities [16, 50]. Ionic impurities, in particular, may be 
expected to have significant effects on the various liquid-liquid, liquid- vapor and liquid-solid 
interfaces; but even in the absence of impurities, an interphase Galvani potential difference 
will appear across the critical liquid-liquid interface together with an associated electrically 
charged double layer [51-54]. 

On the other hand, Moldover [36] in addressing the discrepancy issue notes, as a private 
communication from I. L. Pegg and C. M. Knobler, that further surface-tension measure- 
ments on IBA:W using a different technique (actually light scattering from thermally excited 
capillary waves [50]) yielded the preliminary result K = 11.6 ± 1.0. This value is almost 
three times smaller than found in the HWK experiments with which, as indicated above, 
it cannot be reconciled. Nevertheless, if this value for K is used with the experimentally 
determined values, £^ expt , discussed above it yields a value for S + oc K(£f) 2 close to that 
found experimentally for the other systems and, hence, in accord with the expected universal 
value stated in (A.l) [36, 37]! On this basis one might (despite [42]) be tempted to disregard 
the HWK data and argue that our analysis of the NWW observations should be seriously 
reconsidered. 

Unfortunately, the preliminary study of IBA:W by Pegg and Knobler has never been 
published [50]; nor, to our knowledge, have any subsequent measurements been made by 
others on this system. Thus, while one can and, perhaps, should attribute the nearly three- 
fold difference in the critical surface tension to unidentified impurities or other agencies 
distinguishing the HWK and NWW samples of isobutyric acid from the later sample 
investigated by Pegg and Knobler (in 1984), the true causes of the discrepancies and their 
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significance remain obscure. We must, instead, hope that more sensitive techniques will be 
developed and employed on the IBA:W and other binary systems in order to more stringently 
test the scaling theories for critical endpoints that have been developed. 



[1] J.S. Rowlinson, B. Widom, Molecular Theory of Capillarity, Oxford University Press (1982), 
chap. 9. 

[2] B. Widom, J. Chem. Phys. 67, 872 (1977). 

[3] F. Ramos-Gomez, B. Widom, Physica A 104, 595 (1980). 

[4] N. Nagarajan, WW. Webb, B. Widom, J. Chem. Phys. 77, 5771 (1982), to be denoted here 
as NWW. 

[5] HW. Diehl, M. Smock, Physica A 281, 268 (2000). 
[6] HW. Diehl, M. Smock, Eur. Phys. J. B 21, 567 (2001). 
[7] M.E. Fisher, P.J. Upton, Phys. Rev. Lett. 65, 2402 (1990). 
[8] M.E. Fisher, P.J. Upton, Phys. Rev. Lett. 65, 3405 (1990). 
[9] B. Widom, J. Chem. Phys. 43, 3892, 3898 (1965). 
[10] B. Widom, in Phase Transitions and Critical Phenomena, C. Domb, M.S. Green (Ed.), Aca- 
demic, London (1972), vol. 2, p. 79. 
[11] M.E. Fisher, S.-Y. Zinn, P.J. Upton, Phys. Rev. B 59, 14533 (1999); ibid. 64, 149901 (E) 
(2001). 

[12] S.-Y. Zinn, M.E. Fisher, Phys. Rev. E 71, 011601 (2005) [cond:mat/0410673] to be denoted I. 
[13] H. Kreuser, D. Woermann, J. Chem. Phys. 97, 7757 (1992). 
[14] T. Mainzer-Althof, D. Woermann, Physica A 234, 623 (1997). 

[15] R.G. Howland, N.-C. Wong, CM. Knobler, J. Chem. Phys. 73, 522 (1980) to be denoted 
HWK. 

[16] S.C. Greer, Phys. Rev. A 14, 1770 (1976). 

[17] I.L. Pegg, M.C. Goh, R.L. Scott, CM. Knobler, Phys. Rev. Lett. 55, 2320 (1985). 
[18] M. Amara, M. Privat, R. Bennes, E. Tronel-Peyroz, Europhys. Lett. 16, 153 (1991). 
[19] J. Ataiyan, D. Woermann, Pure Appl. Chem. 67, 889 (1995). 
[20] M.E. Fisher, S.-Y. Zinn, J. Phys. A 31, L629 (1998). 



27 



[21] P. Butera, M. Comi, Phys. Rev. B 62, 14837 (2000). 
[22] A. Pelissetto, E. Vicari, Phys. Repts. 368, 549 (2002). 
[23] R. Guida, J. Zinn-Justin, J. Phys. A 31, 8103 (1998). 

[24] M. Campostrini, A. Pelissetto, P. Rossi, E. Vicari, Phys. Rev. E 65, 066127 (2002). 

[25] S.-Y. Zinn, Ph.D. thesis, University of Maryland, 1997, chap. 5. Note that the value (3 = 0.328 25 
was adopted here for some of the analysis: that corresponds to the value (3 = 0.328 ±0.004 that 
Greer [16] reports as a good fit to the volume fraction difference on neglecting any correction 
terms to the leading power law. However, within the overall accuracy obtainable, this deviation 
from the value in Table I has no serious consequences here. 

[26] S.-Y. Zinn, M.E. Fisher, Physica A 226, 168 (1996). 

[27] M.E. Fisher, G. Orkoulas, Phys. Rev. Lett. 85, 696 (2000). 

[28] Y.C. Kim, M.E. Fisher, G. Orkoulas, Phys. Rev. E 67, 061506 (2003): see especially Eq. (3.20) 

and prior material. 
[29] B. Widom, J.S. Rowlinson, J. Chem. Phys. 52, 1670 (1970). 
[30] N.D. Mermin, J.J. Rehr, Phys. Rev. Lett. 26, 1155 (1971). 
[31] M.E. Fisher, Phys. Rev. Lett. 34, 1634 (1975). 

[32] In contrast to the NWW data, the coexistence curve measured by Greer [16] is somewhat less 
symmetric: on the basis of her data, the diameter can be represented by ~p(T) ~ 0.993 — 0.066t. 
[33] See also Y.C. Kim, M.E. Fisher, J. Phys. Chem. 105, 11785 (2001). 

[34] NWW supposed the volume fraction 4> w of water in the mixture to be related to the mass 
fraction m w via (f) w = m w /[m w + (1 — m w )/c] where c = 0.970 is a constant that represents the 
ratio of the density of pure isobutyric acid to that of pure water at 20 °C. (B. Widom, private 
communication.). 

[35] Moldover in his review of experiments on many systems [36], quotes K = 31.0 from HWK; 
but as we discuss further in the Appendix, Moldover also reports the value K = 11.6 on the 
basis of a private communication from I. L. Pegg and C. M. Knobler: unfortunately, these two 
values are drastically inconsistent. However, since the first accords well with our analysis of 
both the HWK and NWW data, we proceed, at this point, with the estimate (3.12). 

[36] M.R. Moldover, Phys. Rev. A 31, 1022 (1985). 

[37] H. Chaar, M.R. Moldover, J.W. Schmidt, J. Chem. Phys. 85, 418 (1986). It should be noted 



28 



that in summarizing the data on other binary systems in Table I the authors accept, with- 
out discussion, the "preliminary result" of Pegg and Knobler for the critical surface tension 
amplitude of isobutyric acid + water as cited in Ref. [36]. See further the discussion in the 
Appendix here. 

[38] W. Fenzl, Europhys. Lett. 24, 557 (1993) has discussed the NWW data and challenged the 
validity of Antonow's rule for their experiments. However, as pointed out in the article by Fenzl 
himself (B. Widom, private communication), the direction in which it is argued that the rule 
is violated actually contradicts thermodynamic interfacial stability. In Fig. 2 the tendency for 
the high-density surface tension data to rise up (relative to the parallelogram) and vice-versa 
for the low-density data, indicates a trend (violating stability) that Fenzl may have noticed. 
Otherwise, Fenzl also argues that, contrary to our approach and the theoretical analysis of 
NWW, the critical behavior they observed corresponds to a very small surface field h\ (and, 
correspondingly, small Mi). In light of the observations just made concerning Antonow's rule 
and the successful fits we achieve here, Fenzl's proposal does not seem sustainable. 

[39] A.N. Campbell, E.M. Kartzmark, S.C. Anand, Y. Cheng, H.P. Dzikowski, S.M. Skrynk, Can. 
J. Chem. 46, 2399 (1968). 

[40] LA. McLure, B. Edwards, J. Chem. Phys. 70, 3999 (1979). 

[41] We are indebted to B. Widom for drawing our attention to Refs. [39] and [40]. 

[42] This point has been reconfirmed by C. M. Knobler who has kindly re-examined the original 
records of the HWK experiments and sees no reasonable sources of larger errors. 

[43] B. Chu, S.P Lee, W. Tscharnuter, Phys. Rev. A 7, 353 (1973). 

[44] B. Chu, F.J. Schoenes, W.P Kao, J. Am. Chem. Soc. 90, 3042 (1968). 

[45] D. Beysens, A. Bourgou, P. Calmettes, Phys. Rev. A 26, 3589 (1982). 

[46] M. Hasenbusch, K. Pinn, Physica A 203, 189 (1994); see also M. Hasenbusch, Int. J. Mod. 

Phys. C 12, 911 (2001). 
[47] A.J. Liu, M.E. Fisher, Physica A 156, 35 (1989). 
[48] M.P Gelfand, M.E. Fisher, Physica A 166, 1 (1990). 

[49] See, e.g., S. Dietrich in Phase Transitions and Critical Phenomena, C. Domb, J.L. Lebowitz 
(Ed.), vol. 14, p. 1. 

[50] Private communications from C. M. Knobler and I. L. Pegg. The records of the Pegg-Knobler 



29 



observations have been generously made available to us: we have confirmed that they differ 
by more than a factor of two from the HWK and NWW measurements. 

[51] See J.-N. Aqua, S. Banerjee, M.E. Fisher, arXiv cond-mat/0410692 (2004) and [52-54]. 

[52] M.J. Sparnaay, The Electrical Double Layer, Pergamon Press, Oxford (1972), pp. 4, 338. 

[53] J.O'M. Bockris, A.K.N. Reddy, Modern Electrochemistry, Plenum Press, New York (1970), 
vol. 2, sec. 7.2.5. 

[54] PB. Warren, J. Chem. Phys. 112, 4683 (2000). 



30 



TABLE I: Adopted values of fluid critical exponents. For convenience and consistency with the 
universal amplitude ratios assembled in [20] we have employed the corresponding exponent val- 
ues (satisfying the scaling relations). More recent studies [21-24] lead to our currently preferred 
exponent estimates shown here in parentheses. However, for the present analysis the differences 
are of negligible consequence. We may supplement the relations (1.2) with A = (35 = [3 + 7 and 

7 = (2 - n)v. 
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FIGURE CAPTIONS 

FIG. 1. The coexistence curve of the isobutyric acid and water mixture. The dots represent 
the experimental data of NWW. The inner curve is the fit (3.1) using the amplitude B = 
0.0310 with do and a\ from (3.3) while the outer curve employs B = 0.0318 in accord with 
(3.4); the critical temperature has been taken as T c = 26.310 °C in accord with the estimate 
of NWW. 

FIG. 2. Comparison of surface-tension measurements along the coexistence curve: the open 
circles with error bars represent the NWW data. As explained in the text, the parallel- 
ograms have been drawn using £g| 7 = K\t\^ with the amplitude (3.12) obtained from the 
HWK data. The vertical sides of the parallelograms indicate the ranges in which consis- 
tency between the two experiments can be achieved in accord with Antonow's rule. The 
dashed curve, is based on the fit (4.2) obtained below for the background S (T, 0), while the 
solid curve embodies the values (4.6). See text in Sec. IV A for details. 

FIG. 3. The fit to the surface tension background in zero-field. See text for details of the 
generation of the data points (black dots). The dashed curve is drawn using the estimates 
(4.2). The solid curve represents the fit (4.6). 

FIG. 4. A fit for the surface tension on the critical isotherm as a function of the ordering 
field. The dots joined by the dashed lines represent the NWW isotherms nearest the critical 
temperature that serve as approximate upper and lower bounds. See text for details. 

FIG. 5. Comparison of scaling plots for the vapor-liquid surface tension near the critical 
endpoint of isobutyric acid + water. The solid curves represent the theoretically predicted 
scaling functions S^ifh) using the ameliorated EdGF theory. The dots are derived from the 
NWW data for T > T c and T < T c that fall within the symmetric range \h\ < 2.4 x 10 20 
used in fitting the surface tension background. The other symbols ( x for T > T c and + for 
T < T c ) are some of the NWW data falling outside the symmetric range of h (and hence 
not covered by the background fit). Note that an offset of AS = 100 has been used for 
clarity in the vertical scales for T < T c . See the text for further details. 

FIG. 6. Illustration of predicted crossings of surface-tension isotherms of isobutyric acid + 
water above T c . The isotherms have been computed using the EdGF scaling functions S^{fh) 
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from I with a background of the form (1.10) suitably matched to the experimental data. The 
dotted, solid, dashed, and dot-dashed curves represent isotherms AT = T — T c = 0, 1.0, 2.0, 
and 3.0 °C, respectively. The narrow horizontal and vertical lines mark the critical values of 
£ and p. The short horizontal lines mark the ordering field values up to which the fitting 
was performed. See the text for further details. 
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